Spin-dependent transport for armchair-edge graphene nanoribbons between ferromagnetic leads 
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We theoretically investigate the spin-dependent transport for the system of an armchair-edge graphene 
nanoribbon (AGNR) between two ferromagnetic (FM) leads with arbitrary polarization directions at low temper- 
atures, where a magnetic insulator is deposited on the AGNR to induce an exchange splitting between spin-up 
and -down carriers. By using the standard nonequilibrium Green's function (NGF) technique, it is demonstrated 
that, the spin-resolved transport property for the system depends sensitively on both the width of AGNR and the 
polarization strength of FM leads. The tunneling magnetoresistance (TMR) around zero bias voltage possesses 
a pronounced plateau structure for system with semiconducting 7-AGNR or metallic 8- AGNR in the absence of 
exchange splitting, but this plateau structure for 8-AGNR system is remarkably broader than that for 7-AGNR 
one. Interestingly, the increase of exchange splitting A suppresses the amplitude of the structure for 7-AGNR 
system. However, the TMR is enhanced much for 8-AGNR system under the bias amplitude comparable to 
splitting strength. Further, the current-induced spin transfer torque (STT) for 7-AGNR system is systematically 
larger than that for 8-AGNR one. The findings here suggest the design of GNR-based spintronic devices by 
using a metallic AGNR, but it is more favorable to fabricate a current-controlled magnetic memory element by 
using a semiconducting AGNR. 

PACS numbers: 73.63.-b, 85.35.-p, 72.25.Hg 



I. INTRODUCTION 



Spintronics utilizes the electron spin degree of freedom to 
carry information in electronic devices and functionalities. Its 
central task is to control electron spin for information storage 
and precessing. ' Graphene, a single layer of carbon, has low 
intrinsic spin-orbit and hyperfine couplings,^ long spin diffu- 
sion lengths (~2ju m).^ These features suggest that graphene 
is a promising candidate for the potential applications in spin- 
tronics. Triggered by these findings, a large number of works 
have been devoted to the spin-dependent transport in graphene 
and graphene-FM heterostructures, as prototypical spintronic 
devices. For example, a graphene spin-valve device has been 
successfully fabricated showing that spin-polarized cuiTents 
can be injected by means of magnetic (cobalt, permalloy, etc.) 
electrodes,"*"^ which shows high spin-polarized injection ef- 
ficiency. Meanwhile, a 10% magnetoresistance (MR) is ob- 
served '" in a graphene contacted by two soft magnetic elec- 
trodes. The further experiments have found that a inserted tun- 
nel film at the graphene/electrode interface favors MR which 
can reach up to ~12%." Additionally, it is expected'^ that 
the spin dependence of the electron energy can be also pro- 
duced by a FM insulator substrate which acts as an effec- 
tive magnetic field parallel to the graphene layer '^''''. The- 
oretically, using the tight-binding model, Brey and Fertig'^ 
found that the MR is rather small due to the weak dependence 
of the graphene conductivity on the electronic parameters of 
FM leads. Using a continuous model by NGF, Ding et al}^ 
have demonstrated that the TMR exhibits a cusp around zero 
bias in the absence of external magnetic field for a similar 
device. Recently, special attentions have been paid to zigzag- 
edge graphene nanoribbons (ZGNRs) spin-valve devices since 
ZGNRs exhibit fascinating phenomenon due to the size con- 
finement and the edge state. '^"'^ For example, some theoret- 



ical studies have predicated that ZGNR-based spin-valve de- 
vices have a very large TMR.^""^^ 

On the other hand, the leads (electrodes) for spin-valve de- 
vices in the previous works have been supposed to be two- 
dimensional (of honeycomb- or square-lattice) type, but more 
realistically they should be modeled as three-dimensional 
semi-infinite slabs. Further, the investigation on the TMR for 
GNR-magne tic junctions is still lacking so far, especially, the 
influence of an exchange field on the spin-dependent trans- 
port. Another important phenomenon in spintronics is the 
reverse effect to TMR, i.e., spin-polarized electrons passing 
from the left FM layer into the right layer, where the mag- 
netization deviates the left by an angle, may exert a torque 
to the right FM lead. This effect is the so-called^''-^^ STT 
which has been extensively studied both theoretic ally ^^"^^ 
and experimentally^^'^"' for conventional magnetic junction 
systems. A demonstration of the spin-transfer phenomenon 
is the current-induced magnetic switching, which has been 
confirmed experimentally in spin-valves^^-^'' and magnetic 
tunnel junctions (MTJs).^^ Thus, current-induced magnetic 
switching provides a powerful new tool for the study of 
spin transport in magnetic nanostructures. In addition, it of- 
fers the intriguing possibility of manipulating high-density 
nonvolatile magnetic-device elements, such as magnetore- 
sistive random access memory (MRAM), without applying 
cumbersome magnetic fields. However, the previous works 
on spin-dependent transport through a graphene-based sys- 
tem are mainly focused on the tunnel current and the TMR 
effect. 10-16.20-25 -pjjg investigation on the cuiTent-induced STT 
in such system is sparsely reported. Although we have ad- 
dressed this issue for a FM/graphene/FM device,^' in this pa- 
per we extend it to a FM/AGNR/FM system. 

In this paper, we present a theoretical investigation on the 
fully spin-dependent transport through an AGNR between two 
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FIG. 1 : Scliematic illustration of the system consisting of an AGNR 
between two FM leads, where a FM insulator layer with a control 
gate is deposited on AGNR. A unit cell of AGNR represented by two 
dashed lines contains n numbers of A and B sublattice sites labeled 
as lA, IS, • • •, nA, nB. Two additional hard walls are imposed on 
both edges at chains j=0 and n+ 1 . The magnetic moment of the two 
leads is aligned at a relative angle Q. 



FM leads with arbitrary polarization direction at low temper- 
atures, where a magnetic insulator is deposited on AGNRs 
to induce an exchange splitting. By using the Keldysh NGF 
method,"*"'^' the density of state (DOS), the linear conduc- 
tance, the differential conductance, and the current for the 
system have been calculated separately for spin-up and -down 
channels, and consequently the TMR dependence on applied 
bias has been calculated. Further, the dependence of STT 
on the bias and the polarization angle of the two leads have 
also been examined. It is demonstrated that, the spin-resolved 
transport property of the system depends sensitively on both 
the widths of AGNRs and the polarization strength of FM 
leads. Around zero bias the TMR versus bias voltage pos- 
sesses a pronounced plateau structure for both semiconduct- 
ing 7-AGNR and metallic 8-AGNR systems without exchange 
splitting. The plateau structure for 8-AGNR system is remark- 
ably broader than that for 7-AGNR one. Interestingly, the in- 
crease of the exchange splitting strength suppresses the ampli- 
tude of this structure for 7-AGNR system. However, the TMR 
is enhanced at bias that ranged from -A to A for 8-AGNR sys- 
tem. Further, the current-induced STT for 7-AGNR system is 
systematically larger than that for 8-AGNR system. 

The rest of the paper is organized as follows. In Sec. II, 
we derive the analytical expressions for the spin-dependent 
DOS, the conductance, the current, and the STT, starting from 
the system Hamiltonian by NGF approach. Some numerical 
examples and the discussions for the results are demonstrated 
in Sec. III. Finally, Sec. IV concludes the paper. 



II. MODEL AND METHOD 

The geometry of the system considered in this paper is 
shown in Fig. 1, where an AGNR contains two sublattices 



denoted by A and B. We use n, the number of A(B)-site atoms 
in a unit cell, to denote AGNR with different width. Ideal (per- 
fect) AGNRs with «=3 j-\ with positive integer j are metallic, 
otherwise are semiconducting.'**''^ From the top down, atoms 
in a unit cell are labeled as lA, IB,- ■ ■, nA, nB. As shown in 
Fig. 1, the hard-wall condition is imposed on both edges at 
chains j=0 and n + l. The definition of GNRs in this work 
is in accord with the previous convention, '** i.e., an AGNR is 
identified by the number of carbon zigzag chains forming the 
width of ribbon, and an AGNR with n carbon chains is named 
as n-AGNR, therefore its width is W-(n - 1) V3fl/2, where 
a-1.42 A is the C-C bond length. 

The system under consideration here is composed of a 
semiconducting 7-AGNR or a metallic 8-AGNR connected to 
two FM leads, which can be considered as a MTJ (or a spin- 
valve) device. A layer of magnetic insulator is deposited on 
the top of AGNR sample to induce an exchange splitting be- 
tween spin-up and -down carriers.'^ The gate allows us to con- 
trol the Fermi level locally, i.e., to create a tunable barrier in 
AGNR. In this way, controlling on both charge and spin car- 
rier concentrations can be achieved. The magnetic moment 
Ml of the left FM lead is assumed to be parallel to the y-axis, 
while the moment Mr of the right FM deviates from the y- 
direction by a relative angle 6. The tunneling current flows 
along the x-axis [see Fig. 1]. We assume that the magnitude 
of exchange splitting can be modulated by the magnetic in- 
sulator. It seems a reasonable assumption that the magnetic 
insulator is made of diff'erent materials from FM leads. 

The total Hamiltonian for the system considered reads 

H ^Hg+Hl + Hr+Ht, (1) 

where He describes the central AGNR region, //l(R) is the 
Hamiltonian for the left (right) FM lead, and Hj for the cou- 
pling between AGNR and leads. Here, spin-orbit interaction 
in graphene or GNRs is neglected because it is too weak and 
is of the order of 3-4 meV,'*^ "'^ it opens up a gap of the order 
of 10"^ meV at the Dirac point. In the tight-binding approxi- 
mation, these partial Hamiltonians can be respectively written 
as following: 

He = ^(eo + sA)aJ^.^^aisi,cr + t ^ (a\.^aBj,o- + H.c), (2) 

where a^. ^ (afii.a-) creates (annihilates) an electron on site / 
with spin cr and sublattice index /3=A(B), cr-l (X) represents 
the spin-up (-down) state of electrons, eo is the on-site energy, 
t (x2.75 eV) the nearest-neighbor hopping energy, and A the 
exchange splitting energy induced by the FM insulator on the 
top of AGNR, s=+l stands for the electron spin parallel (-I-) 
or antiparallel (-) to the exchange field. The sum over (/, j) is 
restricted to the nearest-neighbor atoms. Hamiltonian (2) can 
also be rewritten in momentum space 

He = ^(eo + sA)a^_q^a/j,q^ + ^[0(q)al,qo-«B,q<^ + H.c.],(3) 

where 4>(q) - -f[2e'«'°/^cos(-^^v) + e^'^'"] is the structural 
factor with wave-vector in the x-direction within the first Bril- 
louin zone (0<|^.i|<;^/(3fl)) and the discretized wave-vector in 
they-direction^j=--i^^ (m=l,2, ■ ■ ■, n). 
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Lcr 



(4) 



and 



Hr = J][efi(k) - o-Mr cos ejc^^Cks^ - Mr sin Oc^^^cy^^iS) 

k,cr 

are respectively for the left and right lead, where £^j^^-€i(k) - 
eV - (tMi with applied bias voltage V and magnetic mo- 
mentum Ml, eL,(/{)(k) is the single-particle dispersion for the 
left (right) lead, c^^^^ (%l(r)o-) creation (annihilation) 
operator of an electron with wavevector k in the left (right) 
lead, (T denotes the opposite spin-polarization with respect 
to cr. Note that here the coupling between the left lead and 
AGNR involves the B sublattice, while it between the right 
lead and AGNR involves the A sublattice (see Fig. 1). How- 
ever, Hamiltonian (5) will be slightly different if the coupling 
at two interfaces involves the same sublattice, but the result is 
unchanged. 

The coupling Hamiltonian 



kLqo- 



(6) 



ksqa 



where (Ti^q) is the coupling matrix between left (right) 
lead and AGNR, N is the number of sites in sublattice A or B. 
It should be pointed out that the annihilation operator a^^qo- in 
Eq. (6) should be replaced by flA.qo- if the coupling between 
right lead and AGNR involves the A sublattice. 

Now we employ the usually defined retarded, ad- 
vanced and correlated Green's function (GF)'"'''*' 
-;6»(f2 - fi)<{flq^(f2),flqv'(fi)))' 
/0(-f2+fi)<{flq,x(f2),flqv'(fi)))' and 
Gqlq^itiJi) = iXaqv'(?i),aqo-(?2)) to the total Hamil- 
tonian in spin space. When the unperturbed GF gq^qai^) 

of the AGNR is obtained, one can obtain G^" „, (a>) and 

' qa,q fl^ ' 

Gq^q,^(cj) from the standard GF technique and the Keldysh 



<^qlq'«(f2,fi) 



(8) 



where r'^oj) = S^-''(^^)/[l-^r(w)2'"'"(w)] with = 
1/A^2:q4;q„(a>) and ^''•^co) = +j/2[r^qq'(w) + ^r«qq'7Jt] 
is the retarded/advanced self-energy, and S*^(w) = 
/[/£(w)r£qq'(w) + /R(w)/?rfiqq'/?'] is the correlated 
self-energy. In the above equations, the matrix 



R 



cos(0/2) -sin(0/2) ' 
sin(6i/2) cos(6i/2) 
fnR)ioJ')=l/[e^"~^"^'^'^""''^ -I- 1] is the Fenni-Dirac distiibution 
function of the left (right) FM lead, the hnewidth function 

r„qq'(a.)-^ r„qq^,(a» j 

r;^qq,(w)=2;r2:k ^'kaq^M'^^'^ " ""kacr) describes the 
influence of the leads. Here, we consider electrons near the 
Fermi level Ef which contribute predominantly to the tunnel- 
ing. In this case one may assume that the coupling matrix is 
independent of energy and set ^qq>=^- 

Furthermore, the isolated retarded GF for an AGNR is de- 
fined as 



[oj - (eo + * A)] - £± + iri+ 



(9) 



with a positive infinitesimal 77'^, the disper- 
sion relation £'±=±|0(q)| and the wave function 

I^^)±=^(I<A)a ± y[^W)B), where the plus (mi- 
nus) sign applies to the upper (lower) n* (n) band, 

l<^>A= V3tT5 2^^v sm(^m')|A„,) and 

l^>i*= V5tTT^"''=i ^''''^'"'sin(^m')|B.'> with the 
number of unit cells A^^ along the x-direction." A direct 
calculation yields expressions for the GF of an AGNR 



I2a CO — (eo + sA) + irf^ 



NAn + 1) 



2n 



^ n+l ' 



m,m'=l 



[oj - (eo + sA) + - f[l + 4cos2^ + 4cos^cos(3^;t/2)] 



(10) 



n+l 



The exchange term sphts the system into two separate spin obtained as 
subsystems. Therefore, the DOS, the linear conductance, and 

the current for spin-up and -down channels can be respectively ^ 

Pm)(<^) = --^^^amii)^^^' ^^^^ 
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/tu)(^) 



n 



XeV/2 
eV/2 



2k 



(/« - /L)^TT(ii)('^)' 



(12) 



(13) 



where V is the appHed bias vohage, and the matrix elements 
in Eqs. (1 1)-(13) are defined as 

Z—i qa.qa^ \G ,A(jS) G ,XcjS)] rs 



and 



x{u) = [G:;M(/;r^/;')]G:;(w)ri 



.(15) 



On the other hand, the STT is the time evolution rate of 
the total spin of left or right FM lead.^^'^^ Here, we ignored 
the interlayer exchange coupling, so the out-of-plane torque is 
zero.'*'* By means of the NGF method, the current-induced in- 
plane STT along jc' -direction in the {x' ,y') coordinate frame 
[see Fig. 1] exerting on the right FM lead can be obtained^^ 

■(- cos 0o-.r + sin Act,)], (16) 
where ctx and cr, are the Pauli matrices. 

III. RESULTS AND DISCUSSIONS 

In what follows, we present some numerical examples of 
PKX), G|(X), Grf_|(|), TMR, and t^' for the system with 7- or 
8-AGNR at low temperatures according to Eqs. (11)-(16). In 
the calculation, the two FM leads are assumed to be made 
of the same materials, i.e., pi-PR-p, where Pl{R)-{^ur)-\ - 
^L{R)i)/i^L(R)^ + TuR)i) is the polarization sti-ength of the 
left (right) FM. Under the wide bandwidth approximation, 
rz,|(|)=rRm)=ro(l ±/?). The TMR is conventionally defined as 
TMR=(//) - Iap)/Ip, where Ip and I^p are the total current for 
the parallel and antiparallel configuration, respectively. The 
coupling Fo between the AGNR and leads without internal 
magnetization is taken as the energy unit. The chemical po- 
tentials for the two FM leads are set as jJip^Ef +0.5 eV with 
the Fermi energy Ef-Q. The on-site energy eo=0, and Fo^t 
because that they are should be in the same order. 



A. Spin-dependent transport 

In Fig. 2, we present the spin-dependent DOS versus en- 
ergy CO (in units of Fq) in parallel configuration (0=0) with ex- 
change splitting A=0.05 and different polarization strength p 
for the system with 7- and 8-AGNR, respectively. In the case 
of bare AGNR (A=F=0), a zero value plateau in DOS for 7- 
AGNR appears symmetrically with respect to the Fermi level 
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FIG. 2: (Color online) DOS for spin-up [(a) and (c)], and -down 
[(b) and (d)] channels as a function of energy oj (in units of Tq) with 
different polarization p, where the parameters are taken as eo=0, 9=0 
and A=0.05. For comparison, the DOS for bare (A=ro=0) 7-AGNR 
(a) and 8-AGNR (c) is shown in (dark) solid lines. 



[see (dark) solid line in Fig. 2(a)] due to the presence of the 
energy gap for semiconducting 7-AGNR, while for 8-AGNR 
a nonzero symmetrical plateau with respect to the Fermi level 
appears [see (dark) solid line in Fig. 2(c)] because of the gap- 
less energy band in metallic 8-AGNR. However, when a mag- 
netic deposition is applied, the positions of plateaus for both 
7-AGNR and 8-AGNR systems shift with the exchange field 
strength A: the spin-up (-down) DOS deviates the origin point 
by the magnitude of A towards the positive (negative) direc- 
tion of the energy. The reason for this phenomenon is that the 
exchange interaction in the magnetic deposition behaves as an 
effective in-plane magnetic field, and acts on the electrons in 
AGNRs, thus leading to the shift of the energy level related 
to the electron spin. For larger a>, pa- shows an oscillation be- 
havior with sharp peaks, which indicates the buildup of Van 
Hove singularities (VHSs) at subband edge. Additionally, we 
notice that as p increases from zero [(red) dashed lines for 
normal leads] to 0.6 [(green) dash-dotted lines] in Fig. 2, both 
7- and 8-AGNR systems are almost immune from p, which 
is different from the bulk graphene tunneling junction. This 
stems from the fact the VHSs induced by the edge effect of 
AGNRs are independent on p. 

Figure 3 shows the spin-dependent linear conductance (in 
units of e^/h) versus exchange splitting A (in units of Fq) with 
different polarization p in parallel alignment at the Fermi en- 
ergy for the system with 7- and 8-AGNR, respectively. It 
is seen that for both 7-AGNR and 8-AGNR systems, G^ in- 
creases [see Figs. 3(a) and 3(c)] while G| decreases [see Figs. 
3(b) and 3(d)] with the increase of polarization p, which is due 
to the conventional spin-valve effect in the magnetic tunneling 
junction. Namely, with increasing p, the proportion of spin-up 
electrons increases and that of spin-down electrons decreases 
in AGNRs. Consequently, it becomes easier for spin-up elec- 
trons to tunnel through the barrier, but harder for spin-down 



5 



0.4 



"a 0.2 



0.0 



0.0- 




(a)7-AGNF{;spin-up 



J. - - - p^.O 

ll. ) • - 

('.1 ... 



J 



p^'O.3 

vp^'O.e 



1, 1. -I > 



(c)8-AGNR,spJn-ljp \ 

I'll ' 



I.I 



0.5 

A/r. 



1.0 



0.4 



0.2 



0.0 



0.2 



0.0 




{b)7-AGNR, spin-down 



(d)8-AGNR,spin-down 



i \ 

w 



0.5 

A/r 



1.0 




0.4 



&0.2 



0.0 



11 




'•r 




'*' 

/ ■! 




/ * 'j 




/ . 7 





^)7-AGNR,spin-down 



I.:- 



•f- 



|d)8-AGNR, spin-down 



' . i' 



0.0 

-1.0 -0.5 0.0 0.5 1.0 -1.1 
eV/r 



-0.5 0.0 0.5 

eV/r 



1.0 



FIG. 3: (Color online) The linear conductance as a function of ex- 
change splitting A with different polarization p for spin-up [(a) and 
(c)], and spin-down [(b) and (d)] channels. The other parameters are 
the same as those of Fig. 2. 
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FIG. 4: (Color online) The I-V curve with different polarization p 
for spin-up [(a) and (c)], and spin-down [(b) and (d)] chaimels. The 
parameters are the same as those of Fig. 2. 



ones due to the presence of the inverse spin-direction for spin- 
down electrons in the tunneUng process, thus we conclude that 
G| increases while decreases with the increase of p. In ad- 
dition, it is found that for the small A (below the full spin 
polarization'*''), the conductance value is zero for 7-AGNR 
system, while nonzero constant for 8-AGNR system. This is 
related to the presence (absence) of energy gap in semicon- 
ducting 7-AGNR (metallic 8-AGNR). For a larger A (above 
the full spin polarization), the conductance for either 7- or 8- 
AGNR system displays an oscillation enhancement behavior 
with shape peaks, which stems from the resonant tunneling 
through different subbands due to the edge effect. '*^ '*^ 

The bias voltage dependence of current (in units of e/ZzFo) 



FIG. 5: (Color online) The differential conductance as a function 
of bias V with different polarization p for spin-up [(a) and (c)], and 
spin-down [(b) and (d)] channels. The other parameters are the same 
as those of Fig. 2. For comparison, the differential conductance for 
7- and 8-AGNR with A=p=0 is also shown in (dark) solid lines. 



for the 7-AGNR and 8-AGNR systems with different polar- 
ization p under the parallel configuration is demonstrated in 
Fig. 4. It is clearly seen that the cuiTent has a step-like struc- 
ture as a function of bias voltage, which results from the con- 
stant DOS at the Fermi energy and resonance due to VHSs. In 
particular, the dependence of the current on the polarization 
p exhibits different behaviors for the 7-AGNR and 8-AGNR 
systems. With increasing p, the current increases in the entire 
bias voltage for the 8-AGNR system, while for the 7-AGNR 
system, the current rises only in the large bias, and almost 
remains zero at small bias voltage. The explanation for this 
phenomenon is as follows; the DOS for the 7-AGNR system 
vanishes near the Fermi level, thus the 7-AGNR system re- 
sembles a insulator-like barrier at the low energy. In this case, 
the electrons will difficultly tunnel through the 7-AGNR sys- 
tem at low bias. 

In Fig. 5, we demonstrate the dependence of differential 
conductance Gd-dlldV (in units of e^//i) on the bias voltage 
for the system with different polarization p under the parallel 
configuration. For comparison, we plot the differential con- 
ductance of the 7- and 8-AGNR systems for A=/:i=0 in (dark) 
solid lines. The differential conductance exhibits the succes- 
sive oscillation peaks corresponding to the resonant tunneling 
through the edge-induced subbands. In addition, we find that 
there exists a respective plateau structure at small bias for the 
7- and 8-AGNR systems. In the range of the plateau, the con- 
ductance approaches to zero for the 7-AGNR system, while 
a finite constant for 8-AGNR system. An interesting charac- 
teristic is that the plateau structure for the 8-AGNR system is 
broader than that for the 7-AGNR system since 8-AGNR has 
a larger interval between the Fermi energy and the lowest sub- 
band [see Fig. 2]. When A is applied, the peak splits into two 
peaks located at the two sides of the original peak for both 
7-AGNR and 8-AGNR systems, and the position of the peaks 
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FIG. 6: (Color online) The dependence of TMR on bias V for (a) 7- 
AGNR, (b) 8-AGNR system with different polarization p and split- 
ting A. 



depends on the magnitude of the exchange splitting. This be- 
havior is due to the exchange splitting of the edge-induced 
subbands. From Fig. 5, G^/j increases while Gd,i decreases 
for both 7- and 8-AGNR systems with increasing p, which is 
a typical spin-valve effect. 

Figure 6 shows the TMR versus the bias voltage V for 7- 
and 8-AGNR systems with different polarization p and ex- 
change splitting A. Firstly, it is noted that the TMR for the 7- 
and 8-AGNR systems has a high value relative to the conven- 
tional magnetic junction; a pronounced plateau for 7- and 8- 
AGNR systems appears at lower bias due to the constant DOS 
near the Fermi energy [see Fig. 2] . These results are quite con- 
sistent with the previous studies. For high bias, the TMR 
exhibits the successive oscillation peaks resided near the edge- 
induced subbands. The amplitude of these peaks can even ap- 
proach to 1 indicating an perfect spin-valve effect. A high 
TMR for the system with an AGNR has also been obtained in 
Refs. [49,50]. It should be pointed out that the result here is 
different from that of Ref. [23] where the TMR for the AG- 
NRs systems is very small in the vicinity of the Fermi energy. 
The plateau structure for 8-AGNR system is systematically 
broader than that corresponding to 7-AGNR system, which 
may be very useful to overcome a well-known shortcoming 
of MRAM, which is the large decrease in TMR with applied 
bias voltage. Interestingly, the increase of A suppresses the 
amplitude of this structure for 7-AGNR system. However, the 
TMR is enhanced within bias range from -A to A for 8-AGNR 
system. This is caused by the exchange splitting influence on 
the DOS as mentioned in Fig. 2. The results are quantita- 
tively different from our previous work for FM/graphene/FM 
system due to the finite size for AGNRs."'-^'^ When the polar- 
ization changes to p=Q.3, as shown by the (blue) dotted lines 
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FIG. 7: (Color online) The dependence of current-induced STT 
for 7-AGNR on (a) relative orientation angle d at the bias voltage 
V=Q.5Tole, (b) bias voltage V at 9=n/3, and (c) exchange splitting 
A at 6=7r/3 and V^=0.5ro/e with different polarization p. The other 
parameters are the same as in Fig. 2. 



and (green) dash-dotted lines in Fig. 6, the TMR for both 
7-AGNR and 8-AGNR systems becomes obviously smaller 
because the leads become less spin-polarized with decreasing 
p.^^ This result is different from Ref. [15], where the TMR is 
nearly independent of the electronic details of the leads. We 
speculate that the structure in the Ref. [15] is very different 
from the structure considered here. Our results suggest that 
it is more favorable to fabricate a graphene-based spin-valve 
device by using a metallic AGNR. 



B. Spin transfer torque 

Figures 7 and 8 show the current-induced STT (t^ ) as func- 
tions of angular 6, bias voltage and exchange splitting strength 
for 7-AGNR and 8-AGNR systems with different polarization 
p. It is seen that versus the angle 6 shows a sine-like be- 
havior for 7-AGNR and 8-AGNR systems, which is in line 
with the previous findings. ^^'^^ A similar result has also been 
obtained for a FM/NM/FM trilayer system discussed in Ref. 
[52], although the transport mechanisms are different. The 
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FIG. 8: (Color online) The dependence of current-induced STT 
for 8-AGNR on (a) relative orientation angle 9 at the bias voltage 
V=0.5ro/e, (b) bias voltage V at 9=n/3, and (c) exchange splitting 
A at 9=n/3 and V'=0.5ro/e with different polarization p. The other 
parameters are the same as in Fig. 2. 



present result can be easily understood because the spin torque 
is proportional to Sr x (Sl x Sr), where Sl and Sr are the 
spin moments of the left and right FMs, respectively. Thus 
we can conclude that vanishes when the relative alignment 
of magnetization of the two FM leads is parallel (0=0) or an- 
tiparallel {0-n). Furthermore, the STT for 7-AGNR system is 
quantitatively larger than that for 8-AGNR system, suggesting 
that a semiconducting AGNR is a better choice for a current- 
controlled magnetic memory element. 

In addition, it is evident that with increasing polarization p, 
the STT is enhanced, as displayed in Fig. 7(a) and Fig. 8(a). 
This is in line with the statement that is proportional to 
the polarization strength of FM.^^ Similar result is obtained 
in the conventional spin-valve device. From Figs. 7(b) and 
8(b), one can notice the step-like feature in the STT as a func- 
tion of bias voltage for both 7- and 8-AGNR systems, which 
is reminiscent of the current versus the bias. The reason is 
that the in-plane is proportional to the difference of spin- 
polarized current (/■'(0)-/*(;r)),2^'^' where /"(e*) = /|(6i)-/|(0) 
is the spin-current densities for the angle 9. In practice, 
(not shown here) for both 7- and 8-AGNR systems vanishes 
identically due to the absence of minority states available for 
tunneling. Thereby, the bias dependence of remains the 



similar feature to that of /*(0). This is also the reason why the 
anomalous bias dependence of the STT has not been predicted 
in this paper according to Ref . [31], where anomalous bias de- 
pendence of the STT only appears in the case of F{Qi)<mn). 
Another reason for this is that the result is obtained here un- 
der the wide bandwidth approximation. From Figs. 7(c) and 
8(c), it can be found that with the exchange splitting A 
displays a oscillation behavior for the 7-AGNR and 8-AGNR 
systems. This result stems from the fact that the exchange 
splitting shifts the VHSs [cf. Figs. 2(a) and (c)], thus leads to 
the edge-induced subband crossing of the bias windows caus- 
ing the tunneling resonance. 



IV. SUMMARY AND CONCLUSION 

In conclusion, we have demonstrated the effects of the 
AGNR width, the polarization strength of FM leads with ar- 
bitrary polarization directions and the exchange splitting on 
the spin-resolved transport properties for both semiconduct- 
ing and metallic AGNRs. In contrast to other related the- 
oretical works, the leads in our system are not supposed to 
be two-dimensional of honeycomb- or square-lattice type, but 
is more realistical modeled as three-dimensional semi-infinite 
FM slabs. By means of Keldysh NGF method, it is found 
that the spin-resolved transport property for the system de- 
pends sensitively on both the AGNR type and the polariza- 
tion strength of FM leads. Around zero bias the TMR ver- 
sus bias voltage possesses a pronounced plateau structure for 
both semiconducting 7-AGNR and metallic 8-AGNR systems 
without exchange splitting. This phenomenon may be very 
useful to overcome a well-known shortcoming of MRAM, 
which is the large decrease in TMR with applied bias volt- 
age. Remarkably, the plateau structure for 8-AGNR system 
is much broader than that for 7-AGNR system. Interestingly, 
the increase of the exchange splitting A suppresses the ampli- 
tude of this structure for 7-AGNR system. However, the TMR 
is enhanced at bias that ranged from -A to A for 8-AGNR 
system. So it may be useful in the design of spin-valve de- 
vice and graphene nanoribbons-based spintronic devices by 
using a metallic AGNR. In addition, the current-induced STT 
for 7-AGNR system is systematically larger that that for 8- 
AGNR system, which may be more favorable to fabricate a 
current-controlled magnetic memory element by using a semi- 
conducting AGNR. As a development of graphene fabrication 
technology, several tens of nanometer wide GNRs can be ob- 
tained by patterning graphene into a narrow ribbon, and the 
large exchange splitting can also be realizable in GNRs, thus 
the model considered in this paper may have potential appli- 
cation in spintronic devices. 
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